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ABSTRACT 


The  effects  of  strain  hardening,  strain  rate  sensitivity,  thermal  softening, 
heat  conduction  and  the  imposed  strain  rate  on  the  shear  localization  process  in 
plane  strain  compression  arc  examined.  ThC'deformation,  stress  and  temperature 
fields  arc  computed  in  an  infinite  solid  which  contains  a  periodic  rectangular 
array  of  inhomogcncitics.  The  inhomogcncitics  give  rise  to  non-uniform 
deformation  fields  which,  under  certain  conditions,  may  localize  in  the  form  of  a 
shear  band.  Boundary  conditions  arc  prescribed  such  that  the  resulting  fields 
possess  periodicity  with  respect  to  the  inhomogencity  distribution.  In  this 
manner,  attention  may  be  confined  to  a  rectangular  region  of  the  solid  which 
surrounds  a  single  inhomogencity.  Full  two-dimensional  analyses  are  performed 
within  the  context  of  a  viscoplasticity  theory  which,  in  the  rate  independent  limit, 
corresponds  to  flow  theory  with  combined  isotropic  and  kinematic  hardening.  Full 
account  is  taken  of  finite  strain  and  rotation  effects,  but  attention  is  confined  to 
quasi-static  loading.  The  initiation  and  propagation  of  shear  bands  is  examined 
for  the  bounding  theories  of  isotropic  and  kinematic  hardening.  The  predicted 
response  depends  significantly  on  the  multi-axial  hardening  characterization  of 
the  solid.  * 


1.  INTRODUCTION 


Shear  bands  frequently  accompany  the  large  plastic  straining  of  materials. 
The  large  localized  strains  in  a  shear  band  often  precipitate  fracture.  When 
shear  bands  do  not  lead  to  fracture,  the  localized  shearing  greatly  affects 
subsequent  plastic  deformation.  Hence,  shear  bands  have  a  dual 
significance:  as  a  precursor  to  fracture  and  as  a  mechanism  of  large  strain 
plastic  deformation.  In  metals  subject  to  high  rates  of  loading,  shear  bands  can 
form  as  the  result  of  a  thcrmomcchanical  instability.  The  same  metal  may 
also  undergo  shear  localization  at  very  low  strain  rates,  where  thermal  effects 
arc  negligible,  with  some  other  mechanism  causing  the  observed  localization. 

The  important  role  of  thermal  softening  in  triggering  shear  band 
formation  at  high  rates  of  deformation  was  recognized  by  Zener  and  Hollomon 
[1].  Experimental  investigations  of  thermal  softening  induced  shear  bands  have 
been  directed  toward  microstructural  features  contributing  to  and  resulting 
from  localization  and  to  measurements  of  the  strain  and  temperature  conditions 
at  which  the  bands  form,  Rogers  [2],  Wright  and  Batra  [3],  Moss  and  Pond  [4], 
Costin  ct  al.  [5]  and  Hartley  [6].  Analytical  and  numerical  studies  have 
generally  focussed  on  conditions  for  instability  in  simple  shearing  deformations. 
Argon  [7],  Clifton  f 8],  Rccht  [9],  Culver  [10],  Backman  and  Finnegan  [11], 
Merzer  [12],  Frcsscngcas  and  Molinari  [13],  Olsen,  Mcscall  and  Azrin  [14],  Wu 
and  Freund  [15]  and  Shawki  [16],  These  theoretical  studies  have  elucidated  the 
relative  roles  of  thermal  softening,  strain  hardening  and  strain  rate  sensitivity 
on  shear  localization  in  circumstances  where  the  onset  of  localization  docs  not 
involve  a  substantial  change  in  stress  state. 


Shear  band  instabilities  arc  also  observed  at  low  rates  where  thermal 
softening  is  not  significant  and,  furthermore,  in  circumstances  where  there  is  no 
evidence  for  material  softening  due  to  progressive  rupture  on  the  microscale, 
Chang  and  Asaro  [17],  Anand  and  Spitzig  [18],  Larsson  ct  al.  [19],  Under 
isothermal  conditions  and  when  microrupturc  effects  arc  absent,  a  key  feature  of 
plastic  material  response  for  shear  localization  is  the  yield  surface  vertex  structure 
implied  by  the  discrete  nature  of  crystallographic  slip,  Rice  [20],  Nccdlcman  and 
Rice  [21],  Asaro  [22],  Tvergaard  ct  al.  [23]  and  Larsson  ct  al.  [19],  In  this  context 
the  significance  of  a  yield  surface  vertex  lies  in  the  reduced  stiffness  to  a  change 
in  loading  path. 

For  rate  dependent  constitutive  relations,  a  flow  potential  surface  vertex 
is  ruled  out  on  quite  general  gounds,  Rice  [24],  However,  Tvergaard  [25]  and 
Hutchinson  and  Tvergaard  [26]  have  shown  that  a  solid  having  a  smooth  yield 
surface  but  with  a  high  curvature  (relative  to  the  isotropic  hardening  surface) 
at  the  current  loading  point  can  give  predictions  of  necking  and  shear  band 
instabilities  in  line  with  those  based  on  a  corner  theory  of  plasticity. 

In  this  investigation,  we  analyze  the  development  of  localization  from  a 
small  internal  material  inhomogcncity  in  a  situation  where  localization  involves 
a  considerable  change  in  loading  path.  The  formulation  incorporates  the 
effects  of  thermal  softening,  strain  hardening,  strain  rate  sensitivity  and  heat 
conduction.  Full  account  is  taken  of  finite  geometry  changes,  but  inertial 
effects  arc  neglected.  We  confine  attention  to  a  range  of  strain  rates,  from  say 

T  1  "i  1 

10  see'  to  10  see  ,  where  material  rate  sensitivity  is  the  main  time  effect 
and  higher  strain  rates  where  inertial  effects  play  an  important  role  arc 
excluded  from  consideration.  The  particular  boundary  value  problem  analyzed 
is  plane  strain  compression  of  a  block  of  material  containing  a  doubly  periodic 
array  of  inhomogcncitics.  The  doublv  periodic  svmmctrv  permits  the  problem 
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for  the  entire  array  to  be  reduced  to  a  boundary  value  problem  for  one  square 
cell.  No  heat  flux  is  permitted  between  cells  so  that  there  is  an  overall 
adiabatic  response.  However,  within  each  cell  heat  conduction  is  accounted 
for. 

Previous  analyses  of  shear  localization  at  internal  inhomogcncitics  have 
been  carried  out  for  rate  independent  solids  deforming  isothcrmally  by 
Abeyaratne  and  Triantafyllidis  [27]  for  plane  strain  and  by  Freund,  Wu  and 
Toulious  [28]  for  anti-plane  shear.  Abeyaratne  and  Triantafyllidis  [27]  confined 
attention  to  quasi-static  deformations  while  in  Freund,  Wu  and  Toulios  [28]  full 
account  was  taken  of  material  inertia. 

Here,  we  present  results  for  both  isotropically  and  kinematically 
hardening  solids.  A  particular  focus  of  our  work  is  the  influence  of  flow 
potential  surface  curvature  on  the  tendency  for  localization  at  various  strain 
rates  and  we  use  a  kinematic  hardening  flow  rule  to  model  a  "rounded  vertex." 
This  use  of  kinematic  hardening  theory  is  intended  for  loading  paths  that  do 
not  involve  extreme  deviations  from  proportional  loading  and  is  quite  distinct 
in  focus  from  the  use  of  kinematic  hardening  as  a  model  for  Bauschingcr 
effects.  Our  results  exhibit  the  transition  from  plastic  flow  property 
dominated  shear  band  localization  at  low  rates  to  thermal  softening  dominated 
localization  at  higher  strain  rates,  with  the  details  of  this  transition  depending 


sensitively  on  the  description  of  the  flow  potential  surfaces. 


2.1  Incremental  Principle  of  Virtual  Work 

The  governing  equations  for  the  finite  element  discretization  of  the 
momentum  balance  equation  arc  based  on  a  Lagrangian  convcctcd  coordinate 
system.  This  approach  has  been  used  extensively  in  a  variety  of  contexts 
involving  finite  deformation  analysis,  c.g.,  Nccdlcman  [29],  Tvergaard  [30],  and 
for  rate  dependent  constitutive  relations  by  Peirce,  Asaro  and  Nccdlcman  [31]. 
Some  essential  features  of  the  formulation  will  be  outlined  here. 

The  reference  configuration  for  the  Lagrangian  formulation  corresponds  to 
the  initial  unstressed  state  with  the  material  points  identified  by  the  Cartesian 
coordinates  x1.  In  the  current  deformed  state,  the  material  points  arc  at  x1  +  u1 
where  u1  arc  the  components  of  the  displacement  vector  on  the  referential  base 
vectors  and  F  =  3(x+u)/3x  is  the  deformation  gradient. 

A  current  state  of  approximate  equilibrium  is  known  at  time  t,  and  the 
incremental  principle  of  virtual  work  is  used  to  compute  the  time  rate  of  change 
of  the  field  variables  for  a  time  increment  At.  Under  plane  strain  conditions 
with  body  forces  neglected,  this  principle  may  be  expressed  as 

|  ii'j  SFjj  dA  =  }T‘  5ui  dS  +  ~  [  {  T‘  5ui  dS  •  J  n']  SFji  dA]  (2-D 

AS  S  A 

where  A  and  S  denote  the  interior  and  the  boundary  of  the  body  in  the  reference 
configuration,  respectively,  and  (')  indicates  the  material  time  derivative.  The 
components  nlJ  arc  the  contravariant  components  of  the  nominal  stress  tensor;  T1 
arc  the  components  of  the  nominal  traction  vector,  given  by  T1  =  nJ1  Vj,  where  v 
is  the  unit  normal  to  the  boundary  in  the  reference  configuration;  and  Fjj  arc  the 


covariant  components  of  the  deformation  gradient;  the  components  of  each  of 
these  tensors  being  on  the  reference  base  vectors.  The  term  in  brackets  on  the 


right  side  of  (2.1)  serves  to  correct  slight  deviations  from  equilibrium  which  result 
from  incremental  time  integration. 


2.2  Balance  of  Energy 

The  principle  of  energy  balance  requires  that  the  time  rate  of  change  of 
energy  (internal  plus  kinetic)  be  equal  to  the  power  input  arising  from  internal 
and  surface  forces  and  heating.  The  following  simplifying  assumptions  arc  made: 
1)  the  rate  of  change  of  internal  energy  may  be  expressed  as  pCp3T/3t  where  p  is 
the  mass  density  and  cp  is  the  specific  heat  at  constant  pressure;  2)  Fourier’s  law 
of  heat  conduction  applies  -  the  heat  flux  is  proportional  to  the  temperature 
gradient,  with  the  constant  of  proportionality  being  the  thermal  conductivity; 
3)  the  fraction  of  plastic  stress  work  which  is  converted  to  heat  is  X,  where  X  is 
in  the  range  of  0.85  to  0.95  for  many  metals  (Taylor  and  Quinnev  [32]).  These 
assumptions  allow  the  referential  form  of  the  energy  balance  equation  to  be 
written  as 


pcpiT 


(kV0T)  +  XT :  d 


where  VQ  represents  the  gradient  with  respect  to  the  coordinates  in  the  reference 
configuration,  and  dp  is  the  plastic  part  of  the  rate  of  deformation  tensor.  The 
rate  of  deformation  tensor  d  is  the  symmetric  part  of  FF~*. 

The  use  of  tire  referential  gradient  in  the  energy  balance  equation  is 
justified  as  follows.  Suppose  that  the  deformation  gradient  can  be  expressed  as 
the  product  of  an  clastic  and  a  plastic  deformation  gradient,  i.c.  F  =  F.F^  (see 

c  p 

Lee  [33]  and  Rice  [34]).  The  clastic  deformation  gradient  represents  the 
stretching  of  the  atomic  lattice  along  with  any  rigid  body  rotation.  The  plastic 
deformation  gradient  includes  the  flow  of  material  through  the  lattice.  Since 
heat  conduction  is  affected  primarily  by  the  distance  between  atoms,  a  more 
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appropriate  representation  of  the  energy  balance  equation  is  written  with  VQ 

replacing  VQ.  However,  if  the  magnitude  of  the  elastic  strains  is  small,  the  clastic 
deformation  gradient  is  approximately  a  rotation.  In  this  ease,  assuming  that  the  heat 
conduction  properties  arc  isotropic,  the  form  of  the  energy  balance  equation  given  by 
(2.2)  is  obtained  when  terms  of  order  strcss/elastic  modulus  arc  neglected. 

The  variational  equation  which  will  be  used  in  the  finite  element 
formulation  is  obtained  by  multiplying  (2.2)  by  6T  and  integrating  over  the 
volume  in  the  reference  configuration,  and  is  expressed  by 

p  cp  ^  5T  dA  =  -  |  k(V0T)  (V06T)dA+  X  j  T  :  dP  5T  dA  +  j  k(VQT  v)6T  dS. 
A  A  AS  (2.3) 


Two  dimensionless  groups  characterize  the  thermal  response. 


X  on 


l  = 


f^-p^n^o 


n  = 


Pc  pTo 


(2.4) 


Here,  € n  is  a  reference  strain  rate,  TQ  is  a  reference  temperature,  hQ  is  a 

reference  length  and  oQ  is  a  reference  stress.  The  parameter  $  is  the  Fourier 

modulus,  and  represents  the  ratio  of  heat  dissipated  by  conduction  to  internal  heat 

u 

storage.  The  quantity  (k/pcpen)  defines  a  length  scale  over  which  heat 
conduction  effects  arc  significant.  The  parameter  n  measures  the  magnitude  of 
internal  heat  generation. 


2.3  Boundary  Conditions 

We  consider  a  block  of  material  containing  a  doubly  periodic  array  of 
material  inhomogcncitics  subject  to  plane  strain  compression.  As  in  Fig.  1, 
the  array  is  imagined  to  be  divided  into  square  cells  of  dimension  2hQ  by  2hQ, 
with  one  inhomogcncity  centered  in  each  ecll.  We  restrict  consideration  to 
deformations  such  that  the  straight  lines  bounding  each  cell  remain  straight 


after  deformation  (this  is  a  stronger  constraint  than  required  by  periodicity). 
Furthermore,  attention  is  confined  to  deformations  that  preserves  the  mirror 
symmetry  of  the  array  so  that  straight  lines  connecting  centers  of  the 
inhomogcncitics  also  remain  straight. 

The  compression  axis  is  taken  to  be  aligned  with  the  x  direction,  so  that 
the  boundary  conditions  for  the  quadrant  analyzed  arc 


u1  =  0  T2  =  0 

at 

x1  = 

0 

(2.5a) 

u2  =  0  T1  =  0 

at 

II 

(N 

X 

0 

(2.5b) 

u2  =  -€n(h0  -  U)  T1  =  0 

at 

= 

^0 

(2.5c) 

^0 

|  T2  dx2  =  0  u1  =  w 

at 

X1  = 

ho 

(2.5d) 

0 

o 

II 

> 

H 

>° 

on 

S 

(2.5c) 

where  U  is  the  displacement  along 

x2  = 

=  h0. 

defined  to 

be 

positive 

for 

compression,  and  en  is  the  imposed 

nominal 

strain 

rate. 

The 

velocity 

of 

material  points  along  x'  =  hQ,  denoted 

by 

w,  is 

a  quantity 

determined  in 

the 

analysis  as  described  in  Section  4.2. 

The  boundary  condition  (2.5c)  requires 

that 

there 

be  no  heat 

flux  across 

the  boundary  of  the  quadrant  analyzed.  Thus,  the  thermal  boundary  condition 
is  an  overall  adiabatic  one,  in  that  while  there  is  a  heat  flux  within  a  cell, 
there  is  no  heat  flux  between  cells.  This  gives  periodicity  of  the  temperature 
field  with  respect  to  the  inclusion  distribution. 
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3.  CONSTITUTIVE  RELATIONS 

The  flow  potential  surfaces  arc  taken  to  be  concentric  cylinders  centered 
about  a  stress  state  a  in  stress  space.  The  radius  of  the  flow  potential  surface 
is  denoted  by  and  is  given  by 


’  (f  s"s“]’ 


where  S  =  S-a.  The  tensor  S  is  the  Kirchhoff  stress  deviator,  given  by 
S  =  T  -  1  / 3 T  :  I,  where  I  is  the  identity  tensor.  Here,  the  kinematic  and 
isotropic  hardening  plastic  flow  rules  arc  based  on  Kirchhoff  stress  rather  than  on 
Cauchy  stress.  As  long  as  clastic  strains  remain  small,  there  is  little  difference 
between  the  two  formulations. 

The  rate  of  deformation  tensor  is  expressed  as  the  sum  of  clastic  and 
plastic  parts  by 

d  =  dc  +•  dP  .  (3.2) 

The  clastic  part  is  given  by 

dc  =  l/1  :  f  (3.3) 

where  f  is  the  Jaumann  rate  of  Kirchhoff  stress  and  L  is  a  fourth  order 
tensor  representing  the  elastic  stiffness  of  the  material.  The  plastic  part  of 
the  rate  of  deformation  tensor  is  given  by 

p  3?  - 

‘J  =  3^—  S  .  (3.4) 


Here,  €  is  specified  by  the  power  law  relation 


€  =  (. 


1-13  t  1  -T  )]  ht ; 


where  « 0  is  a  reference  strain  rate,  T0  is  a  reference  temperature,  m  is  the 
rate  hardening  exponent,  and  B  specifies  the  thermal  softening 
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charactcristics  of  the  material.  This  assumed  linear  temperature  dependence 
is  representative  of  carbon  steels  in  the  range  0  -  500 °C,  c.g.  Stakcr  [35], 
although  at  elevated  temperatures  the  flow  stress  is  generally  a  non-linear 
function  of  temperature. 

We  employ  combined  isotropic  and  kinematic  hardening,  Gocl  and  Malvern 
[36],  with  h(?)  and  the  evolution  equation  for  a  given  by 


h(7)  =  X  oQ[l  +  €/€0]N  +  (1-X) 


A  ^  f 

&  =  y  (l  -  X)b  d 


where 


f?1m  °o  _  N-l 

b  =  —  [1-(B(T-T0)1  N  —  [I  +  €/e0]  .  (3.8) 

'-6o-1  €o 

Here,  oQ  is  a  reference  stress,  tQ  is  a  reference  strain,  N  is  the  strain 
hardening  exponent  and  X  is  a  constant  ranging  between  zero  and  unity.  Purely 
isotropic  hardening  corresponds  to  X  »  1,  while  purely  kinematic  hardening  is 
given  by  X  =  0.  The  parameter  b,  given  by  (3.8),  is  taken  so  that  at  constant 
strain  rate  and  temperature,  the  response  for  proportional  loading  is 
independent  of  whether  the  change  in  flow  potential  surface  is  described  by 
isotropic  hardening  or  by  kinematic  hardening,  see  Nccdlcman  and  Tvergaard 


An  expression  for  the  Jaumann  rate  of  Kirchhoff  stress  is  obtained  in 
terms  of  rate  quantities  by  combining  (3.2)  through  (3.4), 

f  =  L  :  d  -  — ^ —  L  :  S  .  (3.9) 

2  oF 

Although  there  is  no  explicit  yielding  in  this  formulation,  for  small  m  there 
is  an  effective  yield  point  in  that  plastic  strain  rates  arc  very  small  when  the 
numerator  is  less  than  the  denominator  in  (3.5). 


4.  FINITE  ELEMENT  IMPLEMENTATION 


4.1  Rate  Tangent  Modulus 


Using  the  expression  for  the  Jaumann  rate  of  Kirchhoff  stress  given  by  (3.9) 
in  the  incremental  principle  of  virtual  work  (2.1)  gives  an  explicit  Euler  time 
integration  method.  This  approach  requires  very  small  time  steps  for  numerical 
stability.  In  order  to  increase  the  stable  stepsize,  Peirce  ct  al.  [38]  employ  a 
formulation  which  results  in  a  tangent  modulus  type  method  for  solving  the 
governing  rate  equations.  The  method  is  explicit  in  that  no  iterations  arc  required. 
In  this  method  the  effective  uniaxial  plastic  strain  rate  7  within  an  increment  is 
expressed  as  a  linear  combination  of  its  rates  at  times  t  and  t  +  At  by 


€  =  (1  -  0)€t  +  0£t+At  (4.1) 

where  0  is  a  parameter  ranging  from  zero  to  unity.  The  plastic  strain  rate  at  time 

t  +  At  is  approximated  by 

- 

-  -  Se"  |  8?|  j_  d7  ■ 

<,+  |SFl,6p  *aflt  etS  tTf  '  <42) 

An  expression  for  cF  is  obtained  by  differentiating  (3.1)  and  using  (3.7)  and  (3.9), 


dF  =  -  L  :  S  :  d  -  [3G  +  (1  -  X)b]e 


where  G  is  the  clastic  shear  modulus  and  b  is  given  by  (3.8). 
Combining  (4.1)  through  (4.3)  and  solving  for  7  yields 


€  =  p  +  -  L  :  S  :  d 

2of 
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whcrc 


P 


8gAt _  T 

D  I'  +m[l  -  0  (T  -  T0)] 


(4.5) 


0At€t 

Q  =  - 

D  m  aF 

6At?t  HG  +  (1  -  *)b 

D  =  1  +  -  -  + 

m  L  oF 

Substitution  of  (4.4)  into  the  expression 
stress  given  by  (3.9)  results  in 


(4.6) 

—  N  — ■  ( 1  +T/€0)N_  1_|  .  (4.7) 

h  €q  -1 

for  the  Jaumann  rate  of  Kirchhoff 


$  =  C  .  d  -  —  L  :  S  (4.8) 

20f 

where 

C  =  L  -  ]"  q  (L  :  S)  (L  :  S)  (4.9) 

The  expression  for  $  given  by  (4.8)  may  be  used  to  compute  the  stress  rates  n'J 
needed  in  the  incremental  principle  of  virtual  work  (2.1).  The  details  of  this 
computation  arc  given  by  Nccdlcman  [39], 


4.2  Enforcement  of  the  Periodic  Boundary  Conditions 

The  mixed  finite  element  -  Raylcigh-Ritz  method  developed  by  Tvergaard 
[40]  is  used  to  enforce  boundary  condition  (2.5d).  Applying  this  method,  the 
finite  element  solution  of  the  momentum  balance  equation  is  decomposed  into 
four  sub-problems,  which  arc  chosen  to  insure  that  the  boundary  conditions  arc 
satisfied.  The  finite  element  discretization  of  the  incremental  principle  of 
virtual  work  (2.1)  is  expressed  as 

K|j  Uj  SUj  =  (cj  +  Sj)  SUj  ,  (4-10) 


L2- 


whcrc  K  is  the  tangent  stiffness  matrix,  c  is  the  right  hand  side  which  arises  from 
the  equilibrium  correction  term  in  (2.1),  and  s  is  the  right  hand  side  generated  by 
the  second  term  on  the  right  side  of  the  expression  for  ^  given  by  (4.8). 

The  nodal  velocities  arc  expressed  as  a  linear  combination  of  the  four 
sub-problems. 


>i  -  \  “m  U-' 


(4.1 1) 


where  am  arc  the  coefficients  to  be  determined,  and  the  velocities  LT111'  arc 
the  solutions  to 


K  U-(m)  =  F-(ni) 
ij  ri 


(4.12) 


where  is  the  right  hand  side  for  the  m-th  sub-problem.  The  vectors  U*- ' 

and  arc  taken  to  be  the  velocities  which  arise  from  the  equilibrium 

correction  and  the  stress  rate  right  hand  sides  respectively,  with  homogeneous 
boundary  conditions.  Therefore,  F^  =  c  and  F^2)  =  s.  In  sub-problem  3  a 
unit  horizontal  velocity  is  prescribed  at  x'  =  hQ,  and  in  sub-problem  4  a  unit 
vertical  velocity  is  prescribed  at  x2  =  hQ,  with  F^  =  F^  =  0. 

Once  the  velocities  associated  with  the  sub-problems  have  been  computed, 
(4.11)  may  be  substituted  into  (4.10)  to  obtain 


^  ^mna 


=  P 

m  rn 


(4.13) 


where 


R  =  K  -  U-(ni)  U-<n> 
I'mn  ij  u  J  i 


Pn  =  (C;  +  Sj)  U;(n) 


(4.14) 


(4.15) 


In  order  that  boundary  condition  (2.5c)  be  satisfied,  the  coefficient  ot^  must  be  the 
prescribed  vertical  velocity  along  x"  =  hQ  (see  2.5c).  Therefore,  the  unknown 
coefficients  ocj  (i  =  1,3)  arc  governed  by 


3 

)  Rmn«m  =  Pn  +  <n  <ho  ‘  U>  Rn4  •  <416> 

m=  1 

The  velocity  w  in  (2.5d)  is  given  by  the  coefficient  a-j. 

4.3  Solution  of  the  Coupled  System 

The  differential  equations  resulting  from  the  finite  element  discretization 
of  the  variational  equations  associated  with  the  momentum  and  energy  balance 
equations  ((2.1)  and  (2.3))  arc  of  the  form 

Kjj[X,U,T,a,T]t  Uj  +  Cjj  [X,U,T,a,T]t  Tj  =  Fj[X,U,T,a,T]t  ,  (4.17) 

Mjj[X]  Tj  +  Dij[X]Tj  =  Qi[X,T,T,dP]  .  (4.18) 

The  equations  may  be  effectively  uncoupled  by  introducing  an  estimate  for 

•  c  . 

the  nodal  temperature  rate  at  time  t  +  0At,  say  Tj,  into  (4.17).  This  is 

accomplished  by  expressing  Tj  as  a  quadratic  function  of  time  based  on  its  values 

at  the  previous  three  points  in  time  at  which  the  solution  has  been  obtained,  and 

•  c 

extrapolating  to  obtain  Tj.  Now  the  nodal  velocities  may  be  computed  from  the 
equation 

Kjj  Uj  -  Fj  -  Cjj  T-  .  (4.19) 

Following  the  solution  of  the  nodal  velocities,  the  displacements,  stresses 
and  internal  variables  arc  updated  using  a  simple  linear  incremental  updating 
scheme.  Next,  the  heat  generation  rate  XT  :  dp  is  computed  for  each  element, 
and  the  nodal  temperatures  arc  obtained  from  (4.18).  The  nodal  temperature 


-14- 


ratc  at  time  t  +  0At  is  then  calculated,  and  the  stepsize  is  reduced  if  the  computed 
and  estimated  temperature  rates  at  that  time  differ  by  more  than  an  allowable 
tolerance. 

When  the  solution  of  the  coupled  system  at  time  t  +  At  has  been  completed, 
this  solution  becomes  the  new  current  state.  The  procedure  is  then  repeated  to 
obtain  the  solution  for  the  next  displacement  increment. 


4.4  Finite  Elements 

The  finite  element  meshes  used  in  the  subsequent  numerical  computations 
consist  of  four-node  quadrilateral  elements,  made  up  of  four  "crossed"  constant 
strain  triangular  sub-elements.  In  presenting  results,  the  quadrilateral  is  regarded 
as  the  basic  element  and  when  reporting  values  of  field  quantities  the  average 
value  over  the  triangles  is  associated  with  the  centroid  of  the  quadrilateral. 

Nagtegaal  ct  al.  [41]  have  shown  that  an  element  of  this  type  can 
accomodate  isochoric  deformations.  This  is  of  significance  since  the  plastic 
strain  rate  is  volume  preserving,  so  that  the  total  deformation  at  large  strains 
is  nearly  isochoric.  Another  reason  for  using  the  crosscd-trianglc  elements  in 
the  present  context  stems  from  the  need  to  resolve  localized  shear  bands.  As 
discussed  by  Tvergaard  ct  al.  [23],  a  mesh  composed  of  crossed  triangles  can 
resolve  narrow  shear  bands  in  four  directions  -  parallel  with  either  the  sides 
or  the  diagonals  of  the  element.  If  the  mesh  is  not  oriented  so  that  one  of 
these  directions  coincides  with  the  critical  orientation  for  shear  bands  then 


localization  can  be  significantly  delayed  or  entirely  suppressed,  [23]. 


5- 


5.  NUMERICAL  RESULTS 


5.1  Material  Properties 


The  material  parameters  used  in  the  calculations  arc  chosen  to  be  representative 
of  a  4340  steel  studied  experimentally  by  Hartley  [6].  These  properties  arc 
specified  by  aQ  =  1250  MPa,  E  =  200,000  MPa,  v  =  0.30,  €Q  =  0.003,  N  =  0.08, 
m  =  0.01,  e0  =  0.001  scc'^p  =  7833  kg/m^,  cp  =  465  J/kg°C,  k  =  54  W/m°C(in 
the  analyses  involving  heat  conduction),  B  =  0.0016  per  °C,  and  T0  =  20 °C. 
The  parameter  x  in  (2.2),  which  represents  the  fraction  of  plastic  work  which  is  locally 
dissipated  as  heat,  is  taken  as  unity  in  order  to  maximize  the  effects  of  thermal 
softening.  Non-dimensional  parameters  which  depend  on  aQ  arc  aQ/ E  =  0.00625  and 
n  =  1.1714  in  (2.4).  The  initial  undeformed  state  is  stress  free  and  at  a  uniform 
temperature  specified  by  TQ. 

Whether  shear  band  development  occurs  under  adiabatic  or  uniform 
temperature  conditions  strongly  depends  on  the  value  of  the  dimensionless 
parameter  £  in  (2.4).  A  priori  knowledge  of  the  thcrmomccluinical  behavior  of  the 
solid  allows  for  the  parameter  hQ  to  be  chosen  so  that  a  varied  thermal  response  is 
obtained  within  the  range  of  prescribed  strain  rates.  Culver  [10]  calculated 
temperature  profiles  in  a  finite  one-dimensional  region  for  uniform  heat  input 
corresponding  to  a  uniform  strain  distribution.  The  results  show  that  the  thermal 
response  is  essentially  uniform  for  5  >  5.  As  5  decreases,  heat  conduction  effects 
diminish,  and  for  $  <  0.02  over  one-half  of  the  deforming  region  is 
adiabatic.  For  the  choice  of  hQ  =  1  mm,  the  response  of  the  solid  varies  from 
essentially  adiabatic  at  «n/€0  =  5*10^  (en  =  500  see'*)  to  having  a  nearly 
uniform  temperature  distribution  at  €n/€0  =  5*10J  (en  =  5  see  ). 

An  inhomogcncity  is  prescribed  by  specifying  a  value  of  the  flow  stress  oQ 
for  the  quadrilateral  element  nearest  the  origin  equal  to  eighty  percent  of  the 


value  of  the  flow  stress  for  the  surrounding  elements.  In  addition,  the  value  of  the 
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thcrmal  softening  parameter  0  in  (3.5)  for  this  element  is  taken  as  0.0017  per  °C 
which  is  slightly  higher  than  that  for  the  other  elements.  Therefore,  this  element 
remains  softer  than  the  surrounding  material  throughout  the  deformation  history. 


5.2  Homogeneous  Deformation 

In  order  to  provide  a  background  for  the  finite  element  results,  the 
behavior  of  the  constitutive  relation  in  homogeneous  plane  strain  compression 
is  analyzed.  Also,  the  homogeneous  solution  is  used  to  obtain  an  estimate  of 
the  critical  strain  at  which  localization  will  occur  so  that  the  mesh  can  be 
properly  oriented  for  shear  bands. 

The  homogeneous  solution  is  obtained  for  an  adiabatic  condition.  The  solid 
is  subject  to  the  boundary  conditions  given  by  (2.5),  with  £n/«0  =  5x10-* 
(en  =  500  see'1).  Figs.  2  to  4  show  stress-strain  curves  where  all  material 
properties  remain  constant  except  for  one  whose  value  is  varied.  All  non-varying 
parameters  arc  assigned  the  values  specified  previously. 

Stress-strain  curves  arc  presented  for  various  values  of  the  strain  hardening 
exponent  N  in  Fig.  3.  Denote  the  strain  corresponding  to  the  point  of  maximum 
stress  by  em.  For  values  of  N  below  approximately  0.20,  £m  increases  with 
increasing  N.  However,  for  values  of  N  greater  than  0.20,  the  trend  reverses. 
Increasing  strain  hardening  leads  to  two  competing  effects.  Stress  levels,  and 
hence  heat  dissipation,  arc  increased,  which  is  destabilizing;  on  the  other  hand 
stiffness  is  increased. 

The  effects  of  the  dimensionless  parameter  0  defined  in  (2.4)  arc  shown  in 
Fig.  4.  This  parameter  controls  the  magnitude  of  the  internal  heat  generation  rate 
during  the  deformation.  Its  effect  on  the  localization  strain  will  depend  to  a  great 


extent  on  the  thermal  softening  characteristics  of  the  material.  For  a  given  value 
of  0,  the  strain  €m  decreases  rapidly  as  n  increases. 


An  interesting  feature  of  the  bchavi  :  of  the  constitutive  model  in  a 
homogeneous  deformation  occurs  when  the  material  softens  thermally,  i.c.  when  0 
is  non-zero.  When  0  =  0,  the  stress-strain  response  is  independent  of  X,  the 
mixed  hardening  parameter  in  (3.6)  and  (3.7).  However,  the  response  of  the  solid 
is  not  independent  of  X  if  it  softens  thermally.  When  0  is  positive,  different 
stress-strain  curves  result  for  different  values  of  X,  with  stresses  always  highest  in 
the  ease  of  kinematic  hardening  (X  =  0).  This  difference  is  caused  by  the 
temperature  dependence  of  the  parameter  b  given  by  (3.8),  which  is  chosen  so  that 
the  isotropic  hardening  and  kinematic  hardening  characterizations  give  identical 
responses  in  proportional  loading  at  constant  strain  rate  and  temperature. 

5.3  Adiabatic  Analysis  with  Isotropic  Hardening 

In  this  calculation  the  response  of  an  isotropically  hardening  solid  under 
adiabatic  conditions  is  determined.  Therefore,  £  =  0  in  (2.4)  and  X  =  1  in  (3.6) 
and  (3.7).  The  solid  is  deformed  at  €n/«0  =  5*10^. 

The  finite  element  mesh  is  composed  of  20  rows  of  29  quadrilateral 
elements,  providing  a  total  of  580  quadrilateral  elements.  Displaced  shapes  of  the 

finite  element  mesh  arc  shown  in  Fig.  5  at  various  levels  of  the  normalized 

•  •  2 
compression,  U/hQ,  w+icrc  U  is  the  displacement  of  nodes  along  the  line  x  =  hQ. 

Although  each  quadrilateral  element  is  comprised  of  four  constant  strain  triangles, 

only  the  boundary  of  the  deformed  quadrilaterals  is  shown  in  plots  of  this  type. 

The  overall  pattern  of  deformation  at  U/hQ  =  0.12,  which  is  slightly  before  the 

maximum  load,  is  essentially  uniform,  as  shown  in  Fig.  5(a)  at  this  deformation 

level.  Only  a  slight  non-uniformity  in  the  deformation  occurs  near  the  inclusion. 

Beyond  this  point  thermal  softening  dominates  over  strain  hardening.  At 


U/hQ  =  0.16  sufficient  strain  accumulations  have  occurred  so  that  a  shear  band 
becomes  disccrnablc  in  the  corresponding  deformed  mesh  plot  in  Fig.  5(b).  As  the 
deformation  proceeds  the  shearing  of  the  elements  in  the  band  occurs  rapidly. 
Beyond  U/h0  =  0.18  the  shear  band  is  fully  developed,  and  is  shown  at 
U/hQ  =  0.20  in  Fig.  5(c).  The  shearing  deformation  is  so  intense  in  some  of  the 
quadrilateral  elements  that  they  have  almost  degenerated  into  triangles. 

Plots  of  load  vs.  normalized  compression  arc  shown  in  Fig.  6.  The  force  per 
unit  reference  area  is 


u 

?  =  h~  \  j2  (ho’x,)dx'  ' 


(5.3) 


0 


The  curve  labelled  I  in  Fig.  6  corresponds  to  this  analysis.  The  load  has  been 
normalized  with  respect  to  the  reference  stress  aQ  in  this  figure.  The  maximum 
value  of  P  occurs  at  approximately  U/hQ  =  0.14,  and  a  shear  band  forms  shortly 
thereafter.  As  the  deformation  proceeds,  the  rapid  softening  occurring  in  the 
shear  band  results  in  the  continual  decrease  of  the  load. 

A  more  quantitative  description  of  the  deformation  pattern  is  shown  by  the 
contours  of  constant  maximum  principle  logarithmic  strain,  6,  in  Fig.  7.  The 


plots  in  this  figure  correspond  to  the  same  normalized  compression  levels  for 
which  the  deformed  meshes  were  presented.  At  U/hQ  =  0.12  in  Fig.  7(a)  strains 
vary  between  0.12  and  0.14  throughout  most  of  the  domain,  with  strains  as  high  as 
0.22  occurring  at  the  junction  of  the  inclusion  and  the  surrounding  material.  The 
pattern  of  straining  changes  considerably  after  the  attainment  of  the  maximum 
load,  as  observed  in  Fig.  7(b)  at  U/hQ  =  0.16.  Here  the  strain  in  the  band  is 
approximately  twice  that  in  the  material  on  either  side  of  it,  which  has  continued 
to  deform  in  a  nearly  uniform  manner.  Also  evident  in  this  figure  is  a  section  of 
a  shear  band  extending  through  the  material  near  the  upper  right  section  of  the 
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mesh.  This  is  part  of  a  shear  band  which  propagates  from  the  inclusion  initially 
located  at  (0,2hQ),  Part  of  the  shear  band  propagating  from  the  inclusion  initially 
located  at  (2h0,2h0)  also  passes  through  the  mesh  in  Fig.  7(b),  oriented  in  a 
direction  opposite  to  the  predominant  band.  The  strain  in  this  band  is  only 
slightly  higher  than  in  the  surrounding  material,  providing  an  indication  as  to  the 
total  distance  of  propagation  of  a  shear  band  originating  from  any  inclusion. 
Beyond  this  level  of  compression,  any  further  increment  of  imposed  deformation 
produces  a  large  increase  in  strain  within  the  band  relative  to  the  surrounding 
material,  as  evident  in  Fig.  7(c)  at  U/hQ  =  0.20. 

The  width  of  the  shear  band  is  essentially  constant  throughout  the 
deformation.  The  discretization  of  the  mesh  sets  a  minimum  width  for  the  shear 
band,  which  is  the  element  size.  In  this  analysis,  the  shear  band  forms  over  two 
adjacent  rows  of  quadrilateral  elements,  with  the  most  intense  deformation 
occuring  in  essentially  one  triangular  sub-element  in  each  of  the  quadrilaterals  in 
the  band. 

The  contours  of  constant  temperature  (in  °C  and  at  the  same  deformation  levels 
as  before)  in  Fig.  8  illustrate  the  rapid  localized  heating  that  occurs.  The 
temperatures  computed  in  this  adiabatic  analysis  represent  only  an  upper  bound  on 
the  true  temperature  distribution.  Prior  to  the  attainment  of  the  load  maximum, 
temperatures  range  between  70  and  110°C  as  shown  in  Fig.  8(a)  at  U/hQ  =  0.12. 
However,  at  U/hQ  =  0.16,  the  temperature  in  the  band  has  increased  to  approximately 
200 °C  while  the  temperature  of  the  surrounding  material  has  risen  to  about  100 °C. 
By  the  time  the  material  is  compressed  to  80  percent  of  its  original  height,  a 
maximum  band  temperature  of  over  350°C  has  occurred,  with  only  a  small 
accompanying  temperature  rise  on  either  side  of  the  band. 
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5.4  Heat  Conduction  with  Isotropic  Hardening 

The  effects  of  heat  conduction  on  the  localization  process  arc  examined  in 
this  finite  element  calculation.  The  material  properties  and  the  imposed  strain  rate 
retain  the  same  values  used  in  the  previous  analysis,  with  the  exception  of 
5  =  2.965x10  in  (2.4).  This  value  of  {  will  be  used  in  all  subsequent 
calculations.  The  20  by  29  finite  element  mesh  used  in  the  adiabatic  analysis  will 
also  be  used  here. 

This  analysis  shows  that  heat  conduction  causes  a  significant  delay  in  the 
initiation  of  shear  band  localization  at  en/€0  =  5x10^  (en  =  500  see'1).  The  true 
stress  of  the  elements  in  the  shear  bands  docs  not  attain  a  maximum  value  until 
U/h0  =  0.12,  as  compared  with  U/hQ  =  0.08  in  the  adiabatic  ease.  A  shear  band 
did  not  become  visible  in  a  deformed  mesh  plot  in  the  adiabatic  analysis  until 
U/h0  =  0.16.  In  the  heat  conducting  solid  a  shear  band  is  not  disccrnablc  until 
approximately  U/hQ  =  0.20. 

Deformed  finite  element  meshes  arc  shown  in  Fig.  9.  Shear  bands  in  the 
early  stages  of  formation  arc  shown  in  Fig.  9(a)  at  U/hQ  =  0.22.  The  shear 
bands  arc  significantly  broader  than  the  ones  in  the  adiabatic  ease,  even  as 
they  begin  to  form.  At  this  deformation  level  significant  shearing  is  occurring 
over  nearly  three  rows  of  quadrilateral  elements.  Heat  conducted  out  of  the 
shear  band  cminating  from  the  inclusion  at  the  origin  results  in  a  stiffening 
within  this  band  relative  to  the  adiabatic  ease.  The  redistribution  of  strain 
associated  with  this  relative  stiffening  activates  a  secondary  shear  band, 
oriented  perpendicular  to  the  primary  one,  as  shown  in  Fig.  9(b)  at 
U/h0  =  0.26.  The  substantially  greater  thermal  softening  in  the  band  in  the 
adiabatic  ease  docs  not  permit  the  formation  of  a  secondary  band.  The 
deformed  mesh  in  Fig.  9(c)  corresponds  to  U/hQ  =  0.30.  The  normalized  width 
of  the  shear  bands  (relative  to  the  thermodynamic  length  scale  (k/pCptn)  )  is 


0.30.  In  this 
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0.726  as  they  form  and  increases  steadily  to  1.307  at  U/hQ  = 
calculation,  a  normalized  width  of  unity  corresponds  to  approximately 
172  microns,  so  that  the  predicted  shear  band  width  is  about  225  microns  at 
U/h0  =  0.30.  In  the  adiabatic  analysis,  the  shear  band  width  remained  nearly 
constant  throughout  the  deformation  history  at  about  100  microns. 

The  maximum  principle  logarithmic  strain  rate  for  the  elements  in  the 
shear  bands  in  the  adiabatic  analysis  increased  monotonically  following  the 
formation  of  the  bands.  However,  this  is  not  the  ease  in  the  present  analysis. 
The  rate  of  change  of  maximum  principle  logarithmic  strain  for  the  elements 
in  the  band  emanating  from  the  origin  increases  to  about  i  =  UxlO-’  see'*  at 
U/hQ  =  0.22,  giving  a  strain  rate  concentration  in  the  band  of  3.4.  From  this 
point  until  U/hQ  =  0.30  the  strain  rate  concentration  decreases,  dropping  to 
about  1.6.  Beyond  U/hQ  =  0.30  it  begins  to  increase  again. 

The  normalized  load  for  the  heat  conducting  solid  is  shown  by  curve  II  in 
Fig.  6.  The  load  peaks  at  U/hQ  =  0.1S  and,  as  in  the  adiabatic  ease,  a  shear  band 
forms  shortly  afterwards.  Beyond  this  point  the  load  decreases  and  then 
increases  again,  varying  between  1.84  oQ  and  1.90  oQ. 

Contours  of  maximum  principle  logarithmic  strain  and  temperature  arc 
presented  in  Figs.  10  and  1  1,  respectively.  They  arc  shown  at  the  deformation 
levels  corresponding  to  the  deformed  meshes  in  Fig.  9. 


5.5  Kinematic  Hardening 

This  calculation  is  identical  to  the  previous  one  with  the  exception  that  the 


solid  hardens  kinematically  (X  =  0  in  (3.6)  and  (3.7)). 


Deformed  finite  element  meshes  arc  presented  in  Fig.  12  at  compression 
levels  of  U/h0  =  0.19,  0.21,  and  0,23.  Figs.  12(a)  and  (b)  show  the  mesh  just  prior 
to  and  just  after  the  formation  of  the  shear  band.  Prior  to  the  onset  of 
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localization,  the  strain  and  temperature  distributions  given  by  the  kinematically 
hardening  model  arc  in  close  agreement  with  the  predictions  of  the  isotropically 
hardening  model.  The  stresses  in  the  kinematically  hardening  solid  arc  about  five 
percent  greater  than  in  the  isotropically  hardening  one.  Beyond  this  point,  the  two 
models  give  very  different  results.  The  shear  band  in  the  kinematically  hardening 
solid  develops  rapidly,  and  is  fully  formed  shortly  after  its  initiation,  as  shown  in 
Fig.  12(c).  The  primary  effect  of  kinematic  hardening  is  to  produce  shear  bands 
at  lower  overall  strain  levels  relative  to  the  isotropically  hardening  model. 
Furthermore,  the  shear  band  is  narrow  and  well-defined,  and  its  width  varies  less 
with  increasing  heat  conduction  than  in  the  isotropic  hardening  ease. 

The  load  for  the  kinematically  hardening  solid  is  shown  by  curve  III  in 
Fig.  6.  The  load  in  this  case  is  consistently  greater  than  for  the  isotropically 
hardening  solid,  as  discussed  in  Section  5.2.  An  abrupt  drop  in  the  load  occurs 
just  after  the  attainment  of  the  load  maximum  at  U/hQ  =  0.19.  In  fact,  the  load 
decreases  much  faster  than  in  the  adiabatic  ease.  This  rapid  softening  of  the 
material  provides  evidence  that  strain  rates  in  the  band  arc  quite  high. 

Contours  of  maximum  principle  logarithmic  strain  and  temperature  arc 
presented  in  Figs.  13  and  14,  respectively,  at  compression  levels  corresponding  to 
the  deformed  meshes  in  Fig.  12.  In  contrast  to  the  isotropically  hardening  solid, 
the  strain  rates  in  the  band  increase  monotonically  throughout  the  deformation. 
The  strain  rates  in  the  band  in  this  ease  arc  approximately  three  times  greater 
than  in  the  isotropically  hardening  solid. 

5.6  Analyses  at  bower  Strain  Rates 

The  calculations  which  follow  examine  the  response  of  the  heat  conducting 
solid  at  lower  imposed  strain  rates,  specifically  e n / e 0  =  5*10^  (en  =  50  see"') 


and  €_/£„  =  5x10-’  (e.  =  5  see'  ).  In  these  calculations,  the  values  of  {  In  (2.4) 
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arc  0.2965  and  2.965,  respectively.  At  each  of  these  strain  rates,  the  response  is 
calculated  for  both  isotropic  and  kinematic  hardening.  In  each  of  the  analyses,  all 
material  properties  arc  assigned  the  values  specified  in  Section  5.1. 

Results  for  the  isotropically  hardening  solid  deformed  at  6n/€Q  =  5xl04 
(«n  =  50  see'*)  arc  shown  in  Figs.  15  through  17.  Deformed  finite  element  meshes 
arc  presented  in  Fig.  15  at  compression  levels  of  U/hQ  =  0.31  and  U/hQ  =  0.40. 
The  mesh  consists  of  20  rows  and  40  columns  of  elements.  Corresponding  contour 
plots  of  e  and  T  arc  shown  in  Figs.  16  and  17,  respectively.  What  appears  to  be 
the  beginning  of  a  rather  broad  shear  band  is  evident  in  Fig.  15(a)  at 
U/hQ  =  0.31,  and  the  strain  distribution  in  Fig.  16(a)  provides  supporting 
evidence  that  localization  is  initiating.  However,  as  the  deformation  continues  and 
the  band  broadens  due  to  heat  flow,  the  maximum  principle  logarithmic  strain 
rates  throughout  the  band  never  exceed  1.5en.  What  occurs  is  a  diffuse  mode  of 
shearing,  as  observed  in  the  plot  of  €  at  U/hQ  =  0.40  in  Fig.  16(b).  The 
temperatures  throughout  the  material  arc  increasing  only  slightly  faster  than  they 
would  in  a  homogeneous  deformation. 

Results  for  the  kinematically  hardening  solid  deformed  at  en/e0  =  5xl04 
arc  shown  in  Fig.  18.  The  mesh  used  in  this  calculation  consists  of  20  rows  and  35 
columns  of  elements.  The  strain  rate  in  elements  which  arc  inclined  at 
approximately  45  degrees  to  the  compression  axis  begins  to  increase  rapidly  at 
U/hQ  =  0.28,  and  a  shear  band  becomes  disccrnablc  in  a  deformed  mesh  plot  at 
U/hQ  =  0.31.  The  shear  band  is  well-developed  at  U/hQ  =  0.33,  as  shown  in  the 
deformed  mesh  plot  at  this  compression  level  in  Fig.  18(a).  Corresponding  contours 
of  €  and  T  arc  shown  in  Figs.  18(b)  and  (c),  respectively. 

When  the  isotropically  hardening  solid  is  deformed  at  en/e0  =  5x10^ 
(«n  =  5  see'*),  heat  conduction  effects  arc  so  predominant  that  at  no  time  during 

O 

the  deformation  docs  the  temperature  vary  by  more  than  0.1  C  throughout  the 
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solid.  The  solid  was  deformed  to  U/hQ  =  0.60,  and  showed  no  signs  of  shear 
banding.  Although  the  strains  near  the  inelusion  arc  as  much  as  ten  percent 
greater  than  in  the  surrounding  material,  the  isotropically  hardening  solid  is  quite 
resistant  to  localization. 

A  shear  band  docs  form  in  the  kinematically  hardening  solid  when  it  is 


deformed  at 

€n/€0  =  5*103. 

The 

band  begins 

to  form  at 

U/hQ  =  0.40,  and  is 

well-defined 

at  U/hQ  =  0.44 

as 

shown  in 

Fig.  19(a). 

The  mesh  in  this 

calculation 

consists  of  20 

rows 

and  55  columns  of 

quadrilaterals.  The 

corresponding  strain  and  temperature  distributions  arc  shown  in  Figs.  19(b)  and 
(c),  respectively.  This  calculation  illustrates  that  a  shear  band  can  form 
although  the  temperature  distribution  is  nearly  uniform.  Unlike  the  previous 
analyses  which  involved  kinematic  hardening,  the  shear  band  here  is  broadened 
by  the  effects  of  heat  conduction,  and  extends  over  nearly  four  rows  of 
quadrilateral  elements. 

The  results  of  the  last  two  calculations  (at  €n/c0  =  5*103)  arc  in  accord 
with  the  results  of  Hutchinson  and  Tvergaard’s  [42]  study  of  shear  band  formation 
in  rate  independent  solids  deforming  in  plane  strain.  The  classical  clastic-plastic 
solid  with  a  smooth  yield  surface  is  highly  resistant  to  localization.  The 
flow  potential  surface  of  the  kinematically  hardening  solid  acts  as  a  "rounded" 
vertex  and  permits  localizaiton  at  physically  achievable  strain  levels  [42], 
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6.  CONCLUDING  REMARKS 

Previous  analytical  and  numerical  studies  of  thermal  effects  on  shear  band 
development  have  focussed  on  the  competition  between  thermal  softening  and 
material  hardening  in  circumstances  where  localization  docs  not  involve  a  change  in 
loading  path.  Such  analyses  apply  to  high  rate  torsion  tests,  Costin  ct  al.  [5], 
Hartley  [6],  Duffy  ct  al.  [43],  Indeed  most,  if  not  all,  detailed  experimental  data  on 
thermally  induced  shear  localization  is  available  from  such  tests.  However,  in  this 
investigation  we  have  confined  attention  to  quasi-static  deformations  and  analyzed  a 
model  thcrmomcchanical  problem  where  shear  localization  is  accompanied  by  a 
substantial  shift  in  stress  state.  Even  when  thermal  softening  is  the  dominant 
destabilizing  mechanism,  the  multi-axial  constitutive  response  of  the  material  plays  a 
major  role  in  the  flow  localization  process. 

Within  the  framework  we  have  used,  there  is  a  continuous  transition  from 
thermal  softening  dominated  localization  at  high  rates  to  plastic  material  property 
dominated  localization  at  low  rates.  A  strong  similarity  of  the  shear  localization 
process  under  adiabatic  and  isothermal  conditions  was  observed  by  Olsen,  Azrin 
and  Tsangarakis  [44]  in  4340  steel  and  they  remarked  that  whatever  mechanism 
plays  a  role  under  essentially  isothermal  conditions  undoubtedly  also  contributes 
to  instability  under  adiabatic  conditions.  They  discussed  void  initiation  and 
growth  as  the  destabilizing  factor  under  isothermal  conditions.  While  the 
weakening  induced  by  progressive  micro-rupture  can  play  an  important  role  in 
triggering  localization  under  isothermal  conditions,  localization  can  also  occur 
under  these  conditions  when  there  is  no  prior  progressive  rupture  on  the 
microscale,  Chang  and  Asaro  [17],  Anand  and  Spitzig  [18]  and  Larsson  ct  al.  [19]. 

We  have  used  kinematic  hardening  to  model  the  path  dependent  hardening 
associated  with  a  sharply  curved  flow  potential  surface.  This  use  of  kinematic 
hardening  is  intended  for  loading  paths  that  do  not  involve  extreme  deviations 


-26- 


from  proportional  loading,  such  as  reverse  yielding,  and  is  quite  distinct  in  focus 
from  the  use  of  kinematic  hardening  as  a  model  for  Bauschingcr  effects. 
Furthermore,  for  the  strain  levels  and  deformation  histories  encountered  here,  the 
use  of  the  Jaumann  rate  in  the  evolution  equation  for  a,  (3.7),  is  not  at  issue. 

Our  results  show  a  significant  dependence  of  the  course  of  shear  band 

development  on  the  multi-axial  constitutive  characterization  of  the  material. 
Specifically,  shear  band  development  is  very  sensitive  to  the  curvature  of  the  flow 
potential  surface;  localization  occurs  sooner  and  more  abruptly  in  a  kinematic 
hardening  solid  than  in  an  isotropic  hardening  solid.  The  effect  of  flow 

potential  surface  curvature  on  shear  band  development  found  in  our  calculations 
is  consistent  with  previous  comparisons  of  the  behavior  of  kinematic  hardening 
and  isotropic  hardening  solids  in  localization  phenomena,  c.g.  Tvergaard  [25]  and 
Hutchinson  and  Tvergaard  [26].  In  this  context,  the  main  effect  of  an  increase  in 
curvature  of  the  flow  potential  surface  is  the  reduction  in  the  material’s  stiffness 
in  response  to  an  abrupt  change  in  loading  path. 

The  inclusion  of  heat  conduction  into  the  formulation  introduces  a  material 
dependent  length  scale  so  that  our  analysis  is  capable  of  predicting  shear  band 

thickness.  We  find  that  with  heat  conduction,  the  more  sharply  rounded  flow 

potential  surface,  as  modelled  by  the  kinematic  hardening  flow  rule,  gives  rise  to 
narrower  shear  bands.  For  the  isotropic  hardening  solid,  with  heat  conduction 
accounted  for,  the  shear  band  width  greatly  increases  with  decreasing  strain  rate, 
compare  Figs.  10  and  16,  as  well  as  becoming  much  less  intense.  Another  effect  is 
seen  by  comparing  Figs.  10  and  13.  For  the  kinematic  hardening  solid.  Fig.  13, 
deformation  concentrates  rather  rapidly  into  a  well  defined  shear  band.  On  the 
other  hand,  for  the  isotropically  hardening  solid,  the  concentration  of  deformation  is 
slow  enough  for  heat  conduction  to  induce  some  stiffening  in  the  shear  band.  This 
leads  to  the  formation  of  the  secondary  set  shear  bands  seen  in  Fig.  10. 
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FIGURE  CAPTIONS 


The  undeformed  periodic  array  of  inclusions.  The  shaded  area  in 
this  figure  represents  the  region  of  the  domain  which  is  modelled 
numerically. 


Normalized  true  stress  vs.  logarithmic  strain  for  a  homogeneous 
deformation  in  plane  strain  compression  with  €n/e0  =  ^x10^ 
(€n  =  500  see'1)  and  various  values  of  the  thermal  softening 
parameter  0  in  (3.5).  The  values  of  all  other  material  properties 
arc  specified  in  Section  5.1. 


Normalized  true  stress  vs.  logarithmic  strain  for  a  homogeneous 
deformation  in  plane  strain  compression  with  £n/£Q  =  ^x10^ 

( r n  =  500  see*1)  and  various  values  of  the  strain  hardening  exponent 
N  in  (3.6).  The  values  of  all  other  material  properties  are  specified 
in  Section  5.1. 

Normalized  true  stress  vs.  logarithmic  strain  for  a  homogeneous 
deformation  in  plane  strain  compression  with  €n/£Q  =  5x10^ 

(«n  =  500  see'1)  and  various  values  of  the  non-dimensional 

parameter  h  in  (2.4)  which  measures  the  magnitude  of  the  rate  of 
internal  heat  generation.  The  values  of  all  other  material  properties 
arc  specified  in  Section  5.1. 


Deformed  finite  element  meshes  for  the  ease  where  4n/€0  =  5x10^ 
(«n  =  500  see'1),  with  isotropic  hardening  and  £  =  0  in  (2.4)  being  the 
adiabatic  limit.  The  values  of  all  other  material  properties  arc  specified 
in  Section  5.1.  (a)  U/h0  =  0.12,  (b)U/hQ  =  0.16  and  (c)  U/hQ  =  0.20. 
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Normalizcd  load  vs.  displacement  for  the  eases  where  (I)  en/£0  =  5x10^ 
and  5=0  with  isotropic  hardening;  (II)  en/e0  =  5x10*  and 
5  =  2.965x10'”  with  isotropic  hardening;  and  (III)  €n/€Q  =  5xl0J  and 
5  =  2.965x10  with  kinematic  hardening.  The  Fourier  modulus  5  is 
defined  in  (2.4),  with  5=0  being  the  adiabatic  limit,  and  the  non-zero 
value  of  5  results  from  the  thermal  conductivity  k  =  54  W/m  °C.  The 
values  of  all  other  material  properties  arc  specified  in  Section  5.1. 


Contours  of  maximum  principle  logarithmic  strain  for  the  ease  where 
€n/^0  =  5X10*  (en  =  500  see"*),  with  isotropic  hardening  and  5  =  0  in 
(2.4)  being  the  adiabatic  limit.  The  values  of  all  other  material 
properties  arc  specified  in  Section  5.1.  (a)  U/hQ  =  0.12,  (b)U/hQ  =  0.16 

and  (c)  U/hQ  =  0.20. 

Contours  of  temperature  (in  °C)  for  the  ease  where  en/€0  =  5x10-* 
( e n  =  500  see'1),  with  isotropic  hardening  and  5  *  0  in  (2.4)  being  the 
adiabatic  limit.  The  values  of  all  other  material  properties  arc  specified 
in  Section  5.1.  (a)  U/hQ  =  0.12,  (b)  U/hQ  =  0.16  and  (c)  U/hQ  =  0.20. 

The  initial  temperature  is  20°C. 


Deformed  finite  element  meshes  for  the  ease  where  €n/c0  =  5x10** 
(«n  =  500  see**),  with  isotropic  hardening  and  5  =  2.965x10*”  in  (2.4). 
The  valuer  of  all  other  material  properties  arc  specified  in  Section  5.1. 
(a)  U/hQ  =  0.22,  (b)  U/hQ  =  0.26  and  (c)  U/hQ  =  0.30. 


Contours  of  maximum  principle  logarithmic  strain  for  the  ease  where 
€n/e0  =  5x10**  («n  =  500  see**),  with  isotropic  hardening  and 

5  =  2.965x10*“  in  (2.4).  The  values  of  all  other  material  properties  arc 
specified  in  Section  5.1.  (a)  U/h0  =  0.22,  (b)  U/h0  =  0.26  and  (c) 

U/h0  = 


[HMi] 
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Contours  of  temperature  (in  °C)  for  the  ease  where  en/£0  =  5x10-* 
1 

(en  =  500  see  ),  with  isotropic  hardening  and  £  =  2.965x10'“  in  (2.4). 
The  values  of  all  other  material  properties  arc  specified  in  Section  5.1. 
(a)  U/hQ  =  0.22,  (b)  U/hQ  =  0.26  and  (c)  U/hQ  =  0.30.  The  initial 
temperature  is  20°C. 


Deformed  finite  element  meshes  for  the  ease  where  €n/€0  =  5x10-* 
1  4  -> 

( € n  =  500  see  ),  with  kinematic  hardening  and  5  =  2.965x10'“  in 

(2.4).  The  values  of  all  other  material  properties  arc  specified  in 

Section  5.1.  (a)  U/hQ  =  0.19,  (b)  U/h0  =  0.21  and  (c)  U/hQ  =  0.23. 


Contours  of  maximum  principle  logarithmic  strain  for  the  ease  where 
«n/€0  =  5xl05  (en  =  500  see"1),  with  kinematic  hardening  and 
5  =  2.965x10"“  in  (2.4).  The  values  of  all  other  material  properties  arc 
specified  in  Section  5.1.  (a)  U/hQ  =  0.19,  (b)  U/hQ  =  0.21  and  (c) 

U/hQ  =  0.23. 


Contours  of  temperature  (in  °C)  for  the  ease  where  €n/e0  =  5x10^ 
( t n  =  500  see  ),  with  kinematic  hardening  and  l  =  2.965x10'“  in 
(2.4).  The  values  of  all  other  material  properties  arc  specified  in 
Section  5.1.  (a)  U/hQ  =  0.19,  (b)  U/hQ  =  0.21  and  (c)  U/hQ  =  0.23.  The 
initial  temperature  is  20  °C. 


Deformed  finite  element  meshes  for  the  ease  where  en/€Q  =  5x10^ 
(«n  =  50  see'1),  with  isotropic  hardening  and  {  =  0.2965  in  (2.4).  The 
values  of  all  other  material  properties  arc  specified  in  Section  5.1.  (a) 
U/ho  =  0.31,(b)  U/hQ 


=  0.40. 


Contours  of  maximum  principle  logarithmic  strain  for  the  ease  where 
€ n/ €q  =  5x10^  (en  =  50  see" *),  with  isotropic  hardening  and  l  =  0.2965 
in  (2.4).  The  values  of  all  other  material  properties  arc  specified  in 
Section  5.1.  (a)  U/hQ  =  0.31,  (b)  U/hQ  =  0.40. 

Contours  of  temperature  (in  C)  for  the  ease  where  £n/e0  =  5x10 
(en  =  50  see'*),  with  isotropic  hardening  and  l  =  0.2965  in  (2.4). 
The  values  of  all  other  material  properties  arc  specified  in  Section 
5.1.  (a)  U/hQ  =  0.31;  (b)  U/h0  =  0.40.  The  initial  temperature  is 
20  °C. 

Results  at  U/h0  -  0.33  for  the  ease  where  en/e0  =  5x10^ 

( £ n  =  50  see'1),  with  kinematic  hardening  and  £  =  0.2965  in  (2.4);  (a) 
deformed  finite  element  mesh,  (b)  contours  of  maximum  principle 
logarithmic  strain,  and  (c)  contours  of  temperature  (in  °C).  The  initial 

O 

temperature  is  20  C.  The  values  of  all  other  material  properties  arc 
specified  in  Section  5.1. 

Results  at  U/hQ  =  0.44  for  the  ease  where  en/€0  =  5*10^ 

( € n  =  5  see'*)  and  (  =  2.965  in  (2.4)  with  kinematic  hardening;  (a) 

Deformed  finite  element  mesh,  (b)  contours  of  maximum  principle 
logarithmic  strain,  and  (c)  contours  of  temperature  (in  °C).  The 

initial  temperature  is  20°C.  The  values  of  all  other  material 
properties  arc  specified  in  Section  5.1. 
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FIGURE  18 


